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Abstract 

A new approach to turbulence modeling in unsteady developed flows has recently 
been introduced [l], based on results of rapid distortion theory. The approach involves 
closing the k-1 equations for the organized unsteady component of the flow by modeling 
local uns teadiness as a rapid distortion of the local structure of the parent turbulent flow, 
in terms of an effective strain parameter oc t g [2j. In this paper, the phase-conditioned 
equations of motion are developed to accomodate a new unsteady dissipation model and 
local effects of the slow-relaxation time scale of the parent flow. The model equations 
are tested against measurements of the response of a fully-developed turbulent pipe 
flow to the superposition of sinusoidal streamwise oscillation. Good agreement is found 
between measurements and predictions over a wide range of frequencies of unsteadiness, 
indicating that this approach may be particularly well suited to modeling of unsteady 
turbulent flows which are perturbations about a well characterized mean. 


1. INTRODUCTION 

Turbulent flows which feature organized unsteady bulk motion are present in numer- 
ous aerodynamic, engineering and biological devices. In recent years, these flows have 
received considerable attention within the canonical geometries of pipes, channels, and 
flat-plate boundary layers, under external forcings such as continuous sinusoidal stream- 
wise ones which are well suited to analytical treatment. The boundary-layer studies of 
Karlsson [3], Cousteix et al. [4], and Brereton et al. [5], the channel- flow research of 
Tardu et al. [6], and the pipe-flow experiments of Tu k Ramaprian [7], Shemer et al. 
[8], Mao k Hanratty [9], and Hwang k Brereton [10] have a covered a wide range of 
the par am eter space accessible to laboratory experiments. A general conclusion from 
these studies is that mean flow profiles of velocity and turbulent stresses are scarcely 
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distinguishable from counterparts measured in steady flow, with the same average ex- 
ternal conditions. The average state of the flow is then practically independent of either 
the amplitude or frequency of forcing. This observation suggests modeling approaches 
which treat forced unsteady flows as perturbations about a parent mean flow, the struc- 
ture of which may be deduced from the large body of information on steady turbulent 
wall-bounded flows. 

Modeling of unsteady turbulent flows has typically been pursued by solving unsteady 
evolution equations incorporating closures based on steady flow information. Examples 
include the Reynolds stress and /c-e calculations of Kebede et al. [11], low-Reynolds num- 
ber K-e modeling of Cotton k Ismael [12], quasi-steady /c-e computations of Mankbadi 
& Mobark [13], /c-w models of Wilcox [14], and the multi-scale models of Wilcox [15] 
and Kim k Chen [16], applied to unsteady flows in [17]. These approaches are prob- 
lematic when near-wall and low Reynolds number functions are employed in steady flow 
forms, since the steady flow relations between near-wall turbulent measures and the 
momentary value of u T no longer apply because of variation in phase of velocities with 
distance from the wall. While attempts have been made to calculate unsteady flows us- 
ing near-wall functions based on a momentary local turbulent Reynolds number, these 
approaches to modeling remain under development. 

Modeling efforts which decouple the mean and unsteady components of the velocity 
field include some preliminary channel-flow models proposed by Acharya k Reynolds 
[18] (quasi-steady perturbations on the mean flow, one-equation k closure) and the 
recent model of Mankbadi k Liu [l], based on results of rapid distortion theory applied 
to homogen eons turbulent flow. In the latter model, a k-€ formulation of the oscillatory 
component of the velocity field was closed by modeling the momentary structure of the 
perturbation field (in this case, the local ratio of oscillatory Reynolds stress to oscillatory 
turbulent kinetic energy) as a function of the effective strain a e fr in the perturbed flow. 
The excellent qualitative (though not quantitative) agreement of predicted turbulent 
kinetic energy profiles with measured data prompted a closer examination of this model. 
The results of a revised rapid distortion turbulence model and comparisons with new, 
high-frequency oscillatory pipe-flow data are described in this paper. 


2. EXPERIMENTS 

Profiles of streamwise velocity and turbulent intensity in fully-developed oscillatory 
pipe flow were measured in a recirculating water facility with a translucent test section 
(57mm in diameter, 160 diameters long, assuring fully-developed flow at downstream 
locations). Flow control was achieved by motoring a profiled sleeve around a longitudinal 
slot milled in the pipe near its downstream end, under the control of a laboratory 
computer. The mean flow Reynolds number (referenced to centerline velocity and pipe 
diameter) was 11,700, with a burst frequency of about 1.7 Hz. Forced oscillation at up 
to 10 Hz could be achieved in this apparatus, with temporal variation of phase-averaged 
velocity always a good representation of a sine wave. In these experiments the amplitude 
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of oscillation at the pipe centerline varied between 19% of the mean centerline velocity 
at 0.25 Hz and 11% at 4.0 Hz. Further details are provided in Hwang & Brereton [10]. 

Measurements of streamwise velocity were made using a laser-Doppler anemometer 
with frequency shifting and a counter, which was interfaced with a laboratory computer 
to allow phase-resolved measurements of the instantaneous velocity of the flow. Since 
highly repeatable periodic motion could be imposed by the flow-control apparatus, a 
phase-averaging procedure was adopted for decomposition of flow variables into mean, 
oscillatory and turbulent components. Statistical convergence in data was assumed 
to have been reached when the fractional tolerance (a measure of differences in (uti) 
over the first and second halves of the data set, normalized by the rms level in {u'u')) 
reached 0.1%, which typically required at least 1000 ensembles. Here ( ) denotes the 
phase average of the argument. 


3. PHASE-CONDITIONED K -l EQUATIONS 


In close proximity to the pipe wall, the radial wall-normal distance R — r may be 
replaced by y and a cartesian representation used. The triple decomposition and ap- 
propriate time and phase averaging operations allow any general dependent variable, 
/(x,t), to be expressed as the summed contribution of the three parts: 

/(x,i)=7(x) + 7(x,t) + /'(x,<). (1) 


which are the mean or time-averaged one, the oscillatory or periodic one, and the 
turbulent component respectively. For fully developed swirl-free flow, the constant- 
property Navier-Stokes equations may be decomposed to yield a streamwise oscillatory 


momentum equation: 


du 

dt 


1 dp 
p dx 


—u'v' + u 

dy 


d 2 u 

dy 2 


( 2 ) 


and an oscillatory turbulent kinetic energy equation: 


dn 

~dt 


du — dU ' — di d (l~ , ~ dH\ 

= -u'v’— - u'v'- u'v'— -— -v'p + v'k - v— 1 - 

dy dy dy dy \p dy J 


(3) 


where both ' ' and ~ are oscillatory operators and e is the oscillatory component of 

the homogeneous dissipation rate. If oscillatory diffusion is assumed to be dominated 
by viscous effects near the wall, rather than turbulent transport or pressure work, and 
second-harmonic turbulence production is a negligible proportion of the fundamental, 


the k equation may be simplified as: 


dK 


-du 


dU d 2 K 


— = —u'v'— u'v '— — h v 2 

dt dy dy dy £ 


— £ 


(4) 
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While $jl and -u'v' may be taken from benchmark data on the steady parent flow 

[17, 18, 26], a k ~ uV closure must be devised to be consistent with results of rapid 
distortion theory, together with an Z model equation which is physically plausible and 
of the correct asymptotic near-wall form. 


4. RAPID DISTORTION THEORY CLOSURE 


Although rapid distortion theory was originally proposed to estimate how turbulence 
is distorted during rapid passage through large-scale straining motions [19], more recent 
studies have shown that this linear theory retains a high level of applicability to slowly 
changing turbulent shear flows [20] — nonlinear processes tend to limit the develop- 
ment which would take place under rapid distortion rather than alter the structure of 
turbulence. In this spirit, Maxey [2] demonstrated how steady turbulent flows in pipes 
and channels with local values of the structural parameter —u'v' pu/u' could be inter- 
preted as locallj -uniform shear flows which rapidly distorted turbulence of an initially 
axisymmetric spectrum tensor at the necessary local strain to produce the same turbu- 
lence structure. Thus results of rapid distortion theory could be used to ascribe local 
values of an effective rapid-distortion strain a tB to steady wall-bounded turbulent shear 
flows j structural parameters such as — u'v' j uV and u'v' fH could then be expressed as 
functions of a tB . 

The asymptotic strain relations of: 

da = dU (5) 

dt dy ’ 

from the definition of strain through rapid distortion at short timescales, and: 


a 



( 6 ) 


which recovers the limit of equilibrium 
eddy) timescale, suggest the generalized 
strains: 

da _ 

dt 


flow with T as a local slow-distortion (large 
local relation which satisfies the asymptotic 


dU a 
dy T ’ 


( 7 ) 


which may also be presumed to model the effective strain a tB in perturbations about 
a steady flow [2]. For the case of streamwise oscillations upon a fully developed parent 
flow, the dominant components of the perturbation and mean shear strains are the 
same. It therefore appears attractive to model the former in terms of the latter. 

The k ~ uV closure is formed by first expressing —u'v'/k = F(a«fl), where is 

evaluated directly from experimental data (*.e. Laufer [21], Lawn [22]). For oscillations 
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which result in small perturbations of a eff about its local parent-flow value, with a eff - 
a eff + S <ff , Fiats) may be linearized as: 

F(a eff ) = F(a eff ) + a t s F'(a tB ) (®) 

leading to the linearized k ~ uV and a eff relations: 

— uV = f-(oe ff ) K + 5 eff F'(a eff ) k and — = -^-jr ■ ( 9a ^) 

Rapid distortion theory characterizations may also be used to expres s compon ents of the 
turbulent kinetic energy in terms of k and a local effective strain. If -u’v'f uV = G(a eff ), 
then small perturbation analysis leads to a relation for the oscillatory component of 
streamwise turbulence intensity in the form: 

-uV = G(a eff ) u'u' + G^a^JuV (10) 

enabling ttV to be evaluated from F, G , and -u'v' from the parent flow, after compu- 
tation of k. Typical characterizations of G^ ff ) in fully developed turbulent pipe flow 
are shown in Fig. 1 (from [2]). 

The oscillatory rate of dissipation e is modeled in the form: 


~ ~ V* 

e = const, k -=- 


( 11 ) 


where £ is a local length scale of dissipation of the parent flow. Thus the phase-averaged 
rate of dissipation (e) is the phase-conditioned turbulent kinetic energy (/c) divided by 
the local small-eddy time scale £/V*. In essence, (e) is treated as proportional to the 
time-dependent turbulent kinetic energy (/c) supplied to the small scales, dissipated over 
time scales of the parent flow; e is then the unsteady part of (e). The assumption that 
the relevant dissipation time scales axe those of the parent flow is justified if small scale 
effects are locally isotropic and out of tune with relatively low wavenumber organized 
forcing. The excellent agreement between high-frequency regions of energy spectra in 
forced and unforced boundary layers [23, 24] supports this view. 

The correct quasi-steady near-wall behavior for the 7 model may be incorporated by 
an additive term which revises the time scale according to the proposal of Reynolds [25] . 
If the dissipation time scale £ is taken as constant multiple of the local mixing length 
£^ = -uV/(at//dy) 2 , then the correct near-wall behavior in (c) may be achieved by 

setting 




C 2 * 3 / 2 \ 

?(dV/dyf) 



( 12 ) 


This approach simplifies the closure as a one-equation /c-£ closure in which the distri- 
bution of the dissipation length scale is taken directly from the mean flow. The model 
follows the asymptote of 7 = 0 at high frequencies ((*c) frozen at its mean value of 
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7c), with insufficient time for vortex stretching to reduce the scale of motions to the 
dissipation range and retain phase memory. It also matches the phase of e (plus oscil- 
latory turbulent diffusion and pressure work) to the phase of k — a more questionable 
assumption which can be assessed when full simulation data for unsteady flows become 
available. 

The model constants scale the dissipation length scale 2 of the parent flow and may 
be deduced from near-wall steady flow-data. 

As y — i *• 0 , (e) = ci c 2 (k) s 

From the parent flow, or equivalently the steady channel-flow simulations of Kim et al. 
[26], cj c<i = 1.4 x 10 -4 . The constant c\ is evaluated from the local ratios of turbulence 
production to dissipation and —u'v'/lz in the parent flow, at the point closest to the 
w all at which the near-wall correction to the dissipation length scale becomes negligible. 
From simulations of steady turbulent channel-flow [26], this correction reaches 4% at 

y+ ~ 12 where the ratio of turbulence production to dissipation is 1.3 and —u'v'/k = 

y \3/2 

0.15. Setting the ratio of production to dissipation equal to cj (/c/ — u'v') / yields c\ — 

0.08, with comparable values of c\ found at adjacent locations. In combination with the 
constraint on the si2e of the corrective term, ci should be evaluated as close to the wall 
as possible; dissipation is strongest in the wall region and is best modeled on parent flow 
data from that region. Also, if the assumption that dissipation takes place over a mean- 
flow time scale is not fully justified, this approach may also be interpreted as a quasi- 
steady model for dissipation, dominated by fine scale motions in the immediate vicinity 
of the wall, which adjust so rapidly to low wavenumber forcing that they effectively 
retain their steady equilibrium values (referenced to the wall variable u T ). Away from 
the wall this is not the case; hence the desirability of a near-wall calibration of model 

coefficients. 


5. SOLUTION PROCEDURE 


The solution procedure required evaluation of F(o:, ff ), a tB , —u'v', and -g- at discrete 
y positions from benchmark data describing the parent flow. A local effective slow- 
distortion timescale T(y) could then be deduced from (6) for each y value. A linear 
algebraic model of <***, as a function of y + , combined with the explicit du+/dy + model 
of Liakopoulos [27], led to a smoothly varying model for T(y + ), plotted in Fig. 2. The 
model reproduces physically plausible features such as a reduced response time as the 
wall is approached and a necessary reduction in the average local size of large eddies. 

By expressing the pressure in (2) as the inviscid core pressure, the closed equation 
set may be expressed as: 


du 

~dt 


dUr 


dt 


d ~ d 2 u 

— u ' v ’ + 1/— 2 ; 
dy oy z 


3 k — —du 

-37 = ~ UV ~ 5 T 

dt ay 


~du 

u'v'— — I- v 


d 2 K „ 


dy T dy 2 


— € 
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da^B _ du . 
dt dy T 

& = -wneu/dy) 2 


( 14 ) 


u f tr e 




In order- to separate variables, u, u'v', k, and a eB were expressed as ue l(Jjt , 
kc™ 1 and a tB e iujt (where ~ denotes a complex first-harmonic amplitude). The closed 
equation set (12) was then combined as a system of first-order ordinary differential 

equations for the pairs of amplitudes of the first harmonic components of and 


& which were in phase and 90° out of phase with the far-field forcing (see Mankbadi 
oy 

Sc Liu [1] for details). 

The structural coupling of the perturbed flow to the parent wall-bounded flow sug- 
gested normalization in wall units (referenced to u T of the parent flow) , which led to 
expressing frequency in the dimensionless form: w + = uu/u^. It has physical signifi- 
cance as the ratio of the near-wall viscous timescale of the parent flow to the timescale 
of forced unsteadiness. The equation system could be solved as a two point boundary 
value problem with u and k equal to zero at the wall and with gradients of zero far 
from the wall, using a relaxation procedure with iteration by Newton’s method. Inte- 
gration was from the wall outward and required very high precision within the Stokes 
layer (y < y/i/fu) to preserve accuracy in the full-field solution. Norris Sc Reynolds [25] 
observed the same sensitivity in their steady k — l computations. 


6. EXPERIMENTAL AND COMPUTATIONAL RESULTS 

Experimental and computational results were compared in the near-wall region of 
turbulent pipe flow for values of w + ranging from 0.01 to 0.16. The lower value corre- 
sponded to slow oscillation at which a quasi-steady response is approached, with small 
variations in the phase of measured quantities in the wall-normal direction. The upper 
value described a rapidly oscillating flow, at which the phase and amplitude of u were in 
good agreement with a (quasi-laminar) Stokes solution. Comparisons of the measured 
and computed amplitudes of u, the amplitude of the oscillatory streamwise velocity, 
are shown in Fig. 3. While the high-frequency results are well matched by computa- 
tions, which describe the characteristic overshoot in u quite well, discrepancies between 
model and experimental results grow with increasing wall-normal distance and decreas- 
ing frequency. The complementary model and experimental results for the phase of u 
(plotted as the phase advance relative to the phase of the centerline oscillatory velocity) 
are shown in Fig. 4. These results also exhibit discrepancies at low frequencies, which 
become more severe as the wall is approached 

The amplitudes and phases of uV obtained from this model are shown in Figs. 5 and 
6. The shapes of the |uV| profiles are in good agreement with measurements, though 
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they appear to slightly underestimate the data at high frequencies. Since the effect 
of imperfect phase-averaging (cycle-to-cycle variation) and wall vibration would be to 
increase the apparent measured amplitude, and since these experimental difficulties are 
more likely to be encountered at higher oscillation frequencies, they may explai n so me 
of the observed discrepancies. The shortcoming of the model in predicting |u'u'| at 
w + = o.Ol for y + > 20 is probably due to the absence of an oscillatory turbulent 
diffusion term, which would be of importance for an adequate description of slowly 
changing k away from the wall. Another possibility is that the assumptions of small 
distortions about a mean no longer apply when the amplitude of u reaches 19% of 
the mean. The model’s shortcomings at low frequencies might also account for the 
discrepancies in prediction of the phase and amplitude of u for w + = 0.01 in Figs. 3 and 
4. 

The phase of v/u' indicates good agreement between the model and data for y + > 10, 
though the same caveats about predictions of the u> + = 0.01 data apply. Closer to the 
wall, where the amplitude of uV was underpredicted at high frequencies, discrepancies 
in phase become significantly larger. It is possible that smoothed near-wall l model 
might produce better results, but at the cost of greater empiricism. In this region at 
high frequencies of unsteadiness, the treatment of turbulence as if in a locally uniform, 
rapidly distorted shear flow may be too simplistic in the immediate vicinity of the wall. 
A careful re-evaluation of F(a eS ), G^) and T in this region might then be required 
to improve model performance. 

The use of a local large-eddy relaxation time scale T(y), in preference to using the 
wall value throughout the flow, made little difference to predicted behavior. This obser- 
vation re inf orces the sensitivity of one-equation models to correct near-wall modeling, 
consistent with the observation for steady one-equation models that there is a narrow 
window for assum ed values of eo f° r which computations remain stable and can match 
external conditions. 

While data describing other phase-conditioned measurements of the Reynolds-stress 
tensor are scarce, the dependence of wall shear stress (to) on w"*” has been measured 
by a number of researchers [7, 9, 10, 29]. The collapse of the phase of t 0 against u> 
is well-established and serves as a good prediction test for turbulence models. The 
predicted phase dependence of r 0 is plotted in Fig. 7 against data reported from a 
number of different experiments and it is clear that the agreement is reasonable over all 
but the low frequencies for which experimental data are available. The asymptotic high- 
frequency approach to the 45° lead demanded by the Stokes solution is satisfied, and 
the transition region between — 0.008 and w"*" — 0,03 is predicted quite well, though 
the model shows a brief overshoot beyond 45° before reaching the high u/ + asymptote. 
It is not clear if this represents a plausible physical effect masked by scatter in the data 
or if it is a shortcoming of the localized model, which an additional diffusion term in 
(7) might improve. The dimensionless form of the wall shear amplitude r 0 is plotted as 
a function of w + in Fig. 8 and again shows reasonable qualitative agreement, though 
the overprediction of the model is pronounced. It is possible that the shortcomings 
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of a localized uniform-shear closure applied close to a solid surface contribute to this 
discrepancy in amplitude, though curiously not in phase. 

In summarizing the model’s performance, it can reproduce a large number of features 
of experimental data without any tuning to the unsteady flow, and exceeds appreciably 
the capabilities of quasi-steady models. There are clearly areas in which improvements 
can be made, and opportunities for more extensive testing against more detailed data 
sets. It se ems likely that the absence of a momentary turbulence diffusion model, the 
simplicity of the strain evolution equation (7), and lack of distinction between wall- 
bounded and free uniform shear in rapid-distortion models, probably limit performance 
and should be areas for future refinement. 

7. CONCLUSIONS 

The modeling of unsteady turbulent flows as rapid distortions of a well-characterized 
parent mean flow represents a new and promising means of predicting turbulence in 
unsteady shear layers. This approach has been incorporated within a preliminary set 
of model equations for the oscillatory field for a wall-bounded unsteady turbulent flow 
which yields results in surprisingly good agreement with experimental data. The only 
empirical constants required for the model are obtained directly from the parent flow 
and so require no tuning. Some shortcomings of the model predictions indicate that 
opportunities for improved performance lie in: *) modeling of a ta near boundaries; **) 
development of a eff evolution equations which match intermediate as well as asymptotic 
states; and m) incorporation of non-linear a tB structural models for flows with more 
extensive distortion about their mean. Also, several component equations of this model 
are proposed to match correct asymptotic forms, without specific regard to their inter- 
mediate behavior, and might benefit from smoothing through the addition of diffusion 

terms. 

This modeling approach may also be applicable to non-equilibrium turbulent shear 
flows in which added distortion is principally in a direction other than that of the 
principal shear of the parent flow. In that case, development of more sophisticated 
rapid-distortion structural models might be pursued to build general turbulence models 
for various distortions about initial flow conditions. 
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Fig. 2. Large eddy relaxation time scale T + as a function of y + . 


12 





Fig. 3. Measurements and predictions of |u|/|uoo| as a function of y + . 
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o , experimental data, or = 0.01 
A, experimental data, w + = 0.04 
□ , experimental data, w + = 0 
V, experimental data, w + 


Fig. 4. Measurements and predictions of (<^ — as a function of y + . 
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Fig. 5. Measurements and predictions of |u't?|/(u r |u 00 |) as a function of y + . 
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Fig. 6. Measurements and predictions of (<t> ) as a function of y + 
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